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COMPARISON  OF  EDGE  DIFFRACTION  THEORIES 


1.  INTRODUCTION 

In  this  report  we  shell  exeat ne  Uflatsev's  Physical  Theory  of  Diffrac¬ 
tion  (PTD)  end  Mltznar's  Increaentel  Length  Diffraction  Coefficient  Method 
(ILDC)  for  the  calculation  of  edge  scatter.  These  two  theories  are  very 
similar  (respectively)  to  Keller's  Geoaetrlcal  Theory  of  Diffraction  and 
Michaell's  Method  of  Equivalent  Currents,  and  aa  shown  by  Knott  [1,2,3], 
these  letter  two  theories  differ  froa  the  foraer  only  in  the  way  In  which 
they  do,  or  do  not,  incorporate  physical  optics  (PO)  effects  explicitly  in 
the  solutions.  In  fact,  the  PTD  and  ILDC  algorithms  are  applied  to  each 
edge  of  e  plate  (or  wedge  face),  and  the  total  scatter  is  then  obtained  by 
taking  the  coherent  sua  of  all  these  edge  contributions  with  the  PO  area 
scattering  Integral.  To  this  recipe  one  auat  also  add  the  effects  of 
standing  wave  interaction  affects,  as  will  be  discussed. 

In  our  comparisons  we  shall  use  the  PTD  as  modified  by  Knott  in  [1]  and 
[2],  and  also  the  ILDC  as  described  by  ftiott  in  [2]  and  [3].  The  PO  effects 
will  be  calculated  according  to  vectorized  versions  of  the  formula's  given 
by  Gordon  in  [4],  and  also  described  by  Knott  in  [1]  and  [2].  Finally,  we 
mention  that  Mlchaell  has  published  a  new  version  of  his  theory  in  a  very 
recent  paper  [5],  but  that  according  to  the  evidence  that  Mlchaell  hiaself 
presents,  the  results  do  not  appear  to  be  in  very  good  agreeaent  with  exper¬ 
imental  values. 

2.  SUMMARY  OF  MAIN  RESULTS 

The  reeults  suaaerlzed  below  in  this  section  will  be  illustrated  by 
graphical  reeults  presented  in  the  next  section,  along  with  additional 
coaaentary  and  referencea. 

(1)  Both  the  ILDC  end  PTD  algorithms  contain  an  irremovable  singularity 
(pole)  at  grazing  incidence  in  the  forward  scattering  direction;  i.e.,  the 
algorithms  blow  up  when  the  incidence  vector  is  parallel  to  the  plate  and 
the  bl-etatlc  scattering  angle  is  180  degrees.  This  is  a  defect  which  all 
of  these  edge  diffraction  theories  inherit  froa  the  original  Soaaerfeld 
"exact”  solution  to  scattering  by  a  seal-infinite  plane  or  wedge  [6].  The 
practical  importance  of  this  defect  is  that  it  prohibits  the  calculation  of 
adge-to-edga  interaction  effects  near  grazing  Incidence,  as  will  be  dis¬ 
cussed  below. 

(2)  More  generally,  closed-fora  hand  calculations  reveal  that  the  out- 
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puts  from  the  1LDC  and  PTD  are  equal  when  the  incidence  and  scatter  vectors 
are  perpendicular  to  the  edge. 

(3)  Both  the  ILDC  and  PTD  are  numerically  unstable  near  normal  Inci¬ 
dence  to  a  wedge  face.  Hence  it  has  been  necessary  to  use  double  precision 
In  both  cases. 

(4)  There  are  some  "apparent*  singularities  in  the  PTD  and  ILDC  (where 
a  singularity  in  one  term  is  cancelled  out  by  a  singularity  in  another 
term).  The  resolution  of  these  apparent  sigularltles  requires  a  careful 
evaluation  of  limits,  and  complicates  the  computer  coding  by  requiring  the 
Inclusion  of  special  cases.  Moreover,  the  coding  for  the  ILDC  formulas  is 
further  complicated  by  the  existence  of  expressions  for  angles  which  are 
possibly  imaginary  in  the  general  blstatlc  case.  (Equations  5-52  and  5-53 
in  Ref.  [2],  or  Equations  10  and  11  in  Ref.  [3]). 

(5)  Although  there  may  be  significant  differences  between  the  PTD  and 
ILDC  results  for  the  scatter  from  a  single  edge,  these  differences  tend  to 
disappear  when  the  total  radar  cross  section  (RCS)  of  a  flat  polygonal 
plate  is  estimated  by  adding  the  effects  of  all  the  edges  to  the  PO  surface 
effect. 

(6)  At  certain  polarizations  one  must  also  take  into  account  the  effects 
of  standing  wave  interactions.  However,  to  this  author's  best  knowledge,  a 
successful  calculation  of  such  interaction  effects  has  been  carried  out  in 
only  one  special  case;  viz.,  that  of  a  rectangle  when  the  line  of  incidence 
is  perpendicular  to  two  of  the  edges  and  parallel  to  the  others.  In  this 
special  case  the  PTD  and  ILDC  theories  give  the  same  results. 

To  summarize,  we  conclude  that  although  the  PTD  is  somewhat  easier  to 
implement  numerically  than  the  ILDC,  there  does  not  appear  to  be,  at  pres¬ 
ent,  any  good  reason  to  prefer  one  of  these  methods  over  the  other:  Both 
methods  have  singularities  at  grazing  Incidence;  neither  method  has  been 
successfully  applied  to  the  calculation  of  standing  wave  effects  at  arbi¬ 
trary  incidence;  and  both  dm t hods  produce  exactly  the  same  results  in  the 
special  case  for  which  the  calculation  can  be  carried  out  (and  experimental 
results  exist.  See  the  discussion  below.) 

3.  GRAPHICAL  RESULTS  AND  ADDITIONAL  COMMENTARY 

In  the  figures  THETA  is  the  off-normal  angle  to  the  face  of  a  plate  or 
wedge,  PHI  is  an  azimuthal  angle  which  equals  zero  when  the  line  of  sight  is 
perpendicular  to  an  edge,  and  HANG  is  the  interior  wedge  angle.  Hence  we 
set  WANG  -  0  whan  we  calculate  the  scatter  from  a  single  flat  plate.  Edges 
are  oriented  so  that  an  edge  is  a  leading  edge  when  THETA  ■  -90  deg. ,  and  a 
trailing  edge  when  THETA  ■  490  deg. 

In  all  the  figures  shown  below  we  only  consider  the  case  of  monostatic 
backseat ter. 

The  E(H)  polarization  case  corresponds  to  the  case  when  the  E(H)  vector 
is  parallel  to  the  plate  or  wedge  face. 
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The  first  8  figures  show  the  scatter  from  a  single  edge  as  calculated  by 
the  ILDC  (heavy  line)  and  PTD  (dotted  line)  algorithms.  The  calculations 
are  carried  out  for  two  values  of  HANG  (0  and  135  degrees)  and  four  values 
of  PHI  (0,  30,  45,  and  80  degrees).  In  Figures  1  and  5  PHI  -  0,  and  the  PTD 
and  ILDC  results  are  exactly  the  same.  The  maximum  RCS  predicted  by  both 
theories  for  monostatic  scatter  from  an  edge  of  length  L  is  given  by  the 
frequency  independent  result 

<w  -  l2/  *  • 

For  flat  plates  this  maximum  is  only  attained  when  the  grazing  angle  is 
zero  and  the  line-of-sight  is  perpendicular  to  an  edge.  Also,  the  E-vector 
must  be  parallel  to  a  leading  edge  or  perpendicular  to  a  trailing  edge. 

This  maximum  is  shown  in  Figure  1  for  the  case  of  E-polarlzation.  The 
corresponding  curve  for  H-polarlzation  can  be  obtained  by  reflecting  the 
curve  in  Figure  1  about  the  axis  THETA  ■  0. 

Figure  9  shows  what  happens  when  the  edge  effects  are  added  to  the  PO 
effect  to  obtain  an  estimate  for  the  total  RCS  of  a  unit  square.  In  these 
calculations  the  effects  of  edge-to-edge  interactions  are  not  taken  into 
account.  The  dotted  line  shows  the  PO  result  and  the  heavy  line  shows  the 
PO  +  ILDC  result.  The  corresponding  PO  +  PTD  results  is  exactly  the  same 
and  is  therefore  not  shown.  The  corresponding  calculated  result  for 
H-polarlzatlon  is  also  the  same;  however  it  is  known  that  there  is  a  real 
physical  difference  between  these  two  polarizations  which  can  only  be  repro¬ 
duced  by  taking  into  account  the  effects  of  edge-to-edge  Interactions.  More 
specifically,  the  flat  tails  at  grazing  Incidence  shown  in  Figure  9  are 
known  to  be  physically  accurate  for  E-polarlzation,  but  the  experimental 
curves  for  H-polarization  tail-off  to  zero. 

Experimental  and  theoretical  curves  for  the  scatter  from  flat  plates 
have  been  published  by  Ross  in  [7].  We  are  referring  specifically  to  Ross's 
figures  3,4,5  which  depict  the  monostatic  RCS  of  square  plates  as  a  function 
of  THETA  for  the  case  when  the  line  of  incidence  is  perpendicular  to  an 
edge.  These  figures  are  not  reproduced  here  because  of  their  small  size 
and  fine  detail;  however,  in  the  remaining  figures  of  this  report  we  have 
used  the  same  height-to-wldth  ratio  as  in  the  figures  of  Ross  to  facilitate 
a  comparison  between  Ross's  results  and  ours.  He  also  note  that  the  E-  and 
H-polarlzation  cases  discussed  here  correspond  to  Ross's  "vertical”  and 
"horizontal"  polarizations. 

Ross  obtained  his  theoretical  curves  by  applying  the  Geometrical  Theory 
of  Diffraction  (GTD)  to  the  two  dimensional  problem  of  diffraction  by  a 
semi-infinite  strip,  and  then  applying  a  "scaling  law”  to  reduce  the  RCS  so 
obtained  (having  the  dimensions  of  length)  to  an  RCS  having  the  dimensions 
of  length-squared.  This  method  is  only  applicable  when  the  line-of-sight  is 
perpendicular  to  one  pair  of  edges.  In  this  case  the  standing  wave  effects 
are  represented  as  interactions  between  the  edges,  and  the  more  general  case 
would  seen  to  require  that  the  GTD,  PTD,  or  ILDC  be  recast  into  a  fora 
involving  vertex  diffraction,  as  is  suggested  by  the  generalized  Fermat 
Principle  which  motivates  the  GTD  [8,  9]. 

He  are  frankly  at  a  loss  to  understand  how  Ross  calculated  the  edge-to- 
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edge  inCerecCions  for  E-polarization,  since  in  this  case  the  classical 
diffraction  coefficients  of  Sommerfeld  vanish  for  a  surface  wave  propagating 
from  one  edge  to  the  opposite  edge.  In  other  words,  according  to  our  under¬ 
standing,  there  should  be  no  edge-to-edge  effects  for  this  polarization. 

At  any  rate,  according  to  Ross's  figures  the  edge-to-edge  effects  are  very 
saall  for  E-polarizatlon,  except  near  grazing  Incidence  where  they  produce 
a  singularity. 

Figures  10-12  show  the  results  of  our  calculation  of  edge-to-edge 
effects  by  the  straightforward  use  of  the  PTD  algorithm  in  the  H-polariza- 
tlon  case.  The  paraaeters  (radar  wavelength  and  side  of  square  in  aeters) 
are  the  saae  as  in  the  corresponding  figures  of  Ross. 

The  curves  labeled  "no  interactions’*  show  the  results  of  adding  the  P0 
scatter  to  the  direct  edge  scatter  (with  no  edge-to-edge  interactions). 

Hence,  these  curves  can  also  be  taken  to  represent  the  total  RCS  in  the  E- 
polarization  case  because  (according  to  our  understanding)  there  are  no 
edge-to-edge  interaction  effects  for  E-polarized  waves  and  because  the 
direct  edge  effects  for  both  polarizations  are  exactly  the  same.  These 
curves  are  in  good  agreement  with  Ross's  “experimental'’  curves  for  E-polar- 
lzation,  and  moreover  do  not  contain  the  singularities  at  grazing  incidence 
that  occur  in  Ross's  calculations. 

The  curves  labeled  “1st  order  interactions”  contain  the  effects  of  rays 
which  proceed  from  the  radar,  hit  an  edge,  propagate  across  the  surface  to 
the  opposite  edge,  and  then  return  to  the  radar.  Higher  order  interactions 
correspond  to  rays  which  propagate  back-and-forth  between  two  opposite  edges 
before  returning  to  the  radar,  but  from  the  figures  it  can  be  seen  that  the 
second  order  interactions  have  almost  no  effect.  The  curves  involving 
interactions  contain  a  singularity  at  grazing  incidence  because  of  the 
fundamental  defect  which  these  edge  theories  inherit  froa  the  Soamerfeld 
solution. 

In  the  PTD  calculation  a  wave  propagated  from  an  edge  to  the  opposite 
edge  is  attenuated  (in  amplitude)  by  a  factor  proportional  to  (L/R) ,  where 
L  is  the  edge  width  and  R  is  the  distance  to  the  opposite  edge.  This  a  far- 
field  result  whose  accuracy  requires  that  this  factor  be  saall,  which  is  not 
the  case  for  squares.  In  fact,  we  have  found  that  the  accuracy  of  the  PTD 
results  in  this  case  can  be  improved  by  certain  modifications  for  which  we 
have  not  theoretical  justification.  Pirst,  we  arbitrarily  Increase  the 
propagation  factor  by  a  factor  of  2.  This  modification  Improves  the 
behavior  at  the  inner  sldelobes.  Second,  we  multiply  the  radiation  ampli¬ 
tude  pattern  by  a  tapering  factor  equal  to  cos^  8  across  the  entire  range 
of  0.  More  specifically,  referring  to  the  diffraction  coefficients  f  and  g 
discussed  by  Knott  in  [1,2,3],  we  multiply  the  g  coefficient  by  this  taper. 
This  second  modification  has  a  small  effect  on  the  inner  sldelobes  and 
causes  the  calculated  radiation  pattern  to  vanish  at  the  tails  in  the 
B-polarlzatlon  case.  (It  also  has  a  very  slight  effect  on  the  E-polariza- 
tlon  calculations.)  However,  the  end  lobes  (near  grazing  incidence)  produced 
by  this  method  are  lower  and  narrower  that  the  end  lobes  seen  in  Ross's 
experimental  curves.  The  results  of  these  modifications  are  shown  in 
Figures  13-15. 
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